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FSR at leading and next-to-leading order 
in the radiative return at meson factories * 
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The impact of final-state radiation (FSR) on the radiative return method 
for the extraction of the e + e~ hadronic cross section is discussed in detail. 
Possible experimental tests of the model dependence of FSR are proposed 
for the 7r + 7r~ hadronic final state. The importance of the ir + ir~j final state 
contribution to the muon anomalous magnetic moment is investigated, and 
a method based on the radiative return is proposed to extract these con- 
tributions from data. 



PACS numbers: 13.40.Em,13.40.Ks,13.66.Bc 



1. Introduction to the radiative return method 

Electron-positron annihilation into hadrons in the low energy region is 
crucial for predictions of the hadronic contributions to a„, the anomalous 
magnetic moment of the muon, and to the running of the electromagnetic 
coupling from its value at low energy up to Mz (For reviews see e.g. [1-7] 
the most recent experimental result for is presented in [S]). Measure- 
ments of the hadronic cross section for electron-positron annihilation were 
traditionally performed by varying the beam energy of the collider. The 
<]?- and B-meson factories allow to use the radiative return to explore the 
whole energy region from threshold up to the energy of the collider. Even 
if the photon radiation from the initial state reduces the cross section by 
a factor 0(a/-7r), this is easily compensated by the enormous luminosity of 
these 'factories'. A number of experimental results based on the radiative 

* Presented by H. Czyz at XXVII International Conference of Theoretical 
Physics, 'Matter To The Deepest', Ustrori, 15-21 September 2003, Poland. 
Work supported in part by EC 5-th Framework Program under contracts HPRN- 
CT-2000-00149, and HPRN-CT-2002-00311 (EURIDICE network) and Polish State 
Committee for Scientific Research (KBN) under contract 2 P03B 017 24. 



(1) 



2 



ustron03'czyz printed on February 1, 2008 



return was already published [9-18] and in the near future one can expect 
much more data covering large variety of hadronic final states. 

The radiative return method ^5] (see also [2U\). relies on the following 
factorisation property of the cross section 



(e e — > hadrons + 7) = H(Q , Sl 7 ) <r(e e — ► hadrons, Q z 



dQ 2 dn i 

(1) 

where Q 2 is the invariant mass of the hadronic system, f2 7 denotes the 
photon polar and azimuthal angles, and the function H(Q 2 , fL,) is given by 
QED lepton-photon interactions, thus known in principle with any required 
precision. The formula (^Q) is valid for a photon emitted from initial state 
leptons (ISR) and what is more important similar factorisation formula 
applies for the emission of an arbitrary number of photons or even lepton 
pairs (23 from initial state leptons. Let's forget for a while about FSR. 
In that case one can, by measuring the Q 2 differential cross section of the 
process e + e~ — > hadrons + photons + (possibly) lepton pairs (called (Trr 
from here on) and knowing function the H(Q 2 , ...), extract the value of 
o(e + e~ — ► hadrons). The ... in H(Q 2 , ...) stand for the phase space variables 
of photons and/or lepton pairs. However, in practice the cross section ctrr 
is measured within a given experimental setup, which corresponds to an 
integral over a complicated phase space of photons and / or lepton pairs and 
thus the use of Monte Carlo event generators for extraction of the a(e + e~ — > 
hadrons, Q 2 ) become indispensable. Such Monte Carlo programs (EVA |191 
122] . PHOKHARA (23121123), were and are bein S developed. The analysis 
presented in this paper is based on the results obtained by means of the 
program PHOKHARA 3.0 j2S] (For further extensive discussions of various 
aspects of the radiative return method not covered by this article see [26- 
31], while for related discussion of the scan method look |31l I32j ). 

Further complication arises as the photons and/or lepton pairs are emit- 
ted also from final state charged hadrons and special care has to be taken 
when using the radiative return method. That problem is discussed exten- 
sively in Sections [21 and 0J while Section |21 is devoted to NLO contributions 
to a M . 



2. Leading order FSR 

Leading order (LO) contributions to the process e + e~ — > 7r + 7r~7 are 
schematically (permutations omitted) shown in Fig^ The interference of 
ISR, which leads to a C-odd (C stands for charge conjugation) configuration 
of 7r + 7r~ pair (Fig. la), with the FSR amplitude from Fig. lb, corresponding 
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to C-even 7r + 7r~ configuration, vanishes if a charge symmetric event selection 
is used. It gives rise, however, to charge asymmetries and charge induced 
forward-backward asymmetries. 
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Fig. 2. Relative contribution of FSR with respect to ISR to the inclusive photon 
spectrum at y/s = 1.02 GeV (without and with cuts (multiplied by a factor 10)) 
(a) and y/s = 10.52 GeV (b). 



The FSR contribution itself can be as big as 20% of the ISR at yfs = 
1.02 GeV if no cuts are applied (FigEk). At s/s = 10.52 GeV however, 
a very energetic photon has to be radiated to produce a pair of charged 
pions with invariant mass around the p resonance, and its momentum has 
to be compensated by the momentum of the charged pions. As a result, 
photon and pions are produced back to back, thus leading to a negligible 
contribution from LO FSR (FigEJa). At DA^NE energy FSR has to be 
controlled through suitable event selection, and an example of a possible 
event selection is shown in Figf5^, (lower curve). 
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Fig. 3. Charge asymmetric distributions of the pions in $ (the angle between normal 
to the production plain and the direction of the initial positron) for different values 
of Q 2 and fixed pion polar angles. 



One can easily reduce the FSR contribution to less than 1% and, rely- 
ing on a MC generator, subtract it from the measured cross section, thus 
being able to apply the procedure of the hadronic cross section extraction 
described in the previous section. The FSR contribution is however model 
dependent and one needs an independent experimental check on the ac- 
curacy of the model used. It can be done just by relaxing the cuts and 
measuring various charge asymmetric distributions. As the actual contri- 
bution of FSR to the radiative return cross section, is of the order of 1%, a 
modest 10% accuracy of the model will lead to an error of 0.1%, sufficient 
for any high precision measurement. Some of the tests of the model used in 
EVA and PHOKHARA for FSR (point-like pions and scalar QED (sQED)), 
proposed in ^S], were already done by KLOE [0J, where it was shown that 
the charge asymmetry 

= N* + (0)-N*-{0) 
{ ' N* + (0) + (6) ' 1 ' 

agrees well with the EVA MC |19j . However additional tests are needed to 
assure the accuracy of the model at the required level and comparisons of 
various charge asymmetric distributions between experimental data and MC 
are indispensable. If only pions four-momenta are measured, as done in the 
KLOE experiment at the moment (see |15U16| ). one arrives at distributions 
as shown in FigEJ With 500 pb _1 , collected till now by KLOE, the 0.1 
nb/bin in the plot corresponds to 2000 events per bin. Thus that kind of 
measurement is feasible and tests can be done with the required precision, 
provided systematic errors are small enough. The nontrivial Q 2 and polar 
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angle dependence of ^-distributions provides profound cross checks of the 
tested model. 

3. Next-to-leading order hadronic contributions to muon a M 




Fig. 4. Hadronic leading order (a) and next-to-leading order (hadrons+photon) (b) 
contributions to the muon anomalous magnetic moment. 

The only reliable method of calculation of the hadronic contribution to 
the muon anomalous magnetic moment ciu, available till now in the low 
energy region, is based on dispersion relations, where the measured e + e~ 
hadronic cross section is convoluted with the known kernel function (for 
definitions and recent reviews see [21 El El)- The LO hadronic contribution to 

enters through the graph(s) shown in Fig@Ji. One class of contributions 
at the next-to-leading order (NLO) is shown in FigEb, where an additional 
photon is attached to charged hadron line(s). For the two pion plus photon 
final state one can try to estimate the contribution using point-like pions 
and sQED or just by inserting a quark loop with photonic corrections. As 
shown in J2E1, one cannot rely on a theoretical model, as the values obtained 
vary very much with the chosen model and do depend very strongly on the 
value of the quark mass. The values obtained are: <5a^ (quark, 7, m q = 
180 MeV) = 1.880 x 10~ 10 , 5a M (quark, 7, m q = 66 MeV) = 8.577 x 10~ 10 
and Sa^TT+ir^ ,j) = 4.309 x 10~ 10 , where m q = 180 MeV describes well 
the LO contributions to , while m q = 66 MeV describes well the lowest 
order contribution to a[Mz)- 

As the values are of the order of the present error of the hadronic con- 
tribution to [5] it is important to measure also the relevant parts of the 
e + e~ — ► 7r + 7r~7 cross section and then, using dispersion relations, get the 
contribution of the 7r + 7r~7 intermediate state to a^. 

This is by no means easy (for more extensive discussions see [221)) as 
the two body cuts of the diagram in Fig@Ja correspond to the radiative 
corrections to the 7T7T7 vertex (FigEJa), while the three body cuts correspond 
to a part of the leading order e + e~ — > tt + tt~j process (FigEt), with the 



G 
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Fig. 5. Leading order (a) and 7r + 7r~7* vertex corrections (b) contributions to the 
reaction e + e~ — ► n + Tr~; leading order FSR contribution to the reaction e + e~ — * 
7r + 7r _ 7 (c). 

photon emitted from the final states only. In practice, when measuring the 
cross section a(e + e~ — ► 7r + 7r~), one always measures in combination the 
contributions from FigOi and |3) plus the 'soft' part of the contribution 
from FigOJ, with the latter depending on the actual experimental setup. 

The virtual plus soft corrections to the tt + tt~j* vertex are negative. As 
a result, the actual contribution from the 'hard' part (photon with energy 
above E cu t) can be larger than the inclusive sum of 'virtual'+'soft'+'hard' 
contributions. This is shown in FigHJ,, where differential contributions to 
a had,7 f rom 7r + 7r -^, intermediate states are compared with the contribution 
from 7r + 7r~. The contribution to a^, integrated over the whole s-spectrum, 




s(GeV 2 ) E cnt [MeV] 



Fig. 6. Differential contribution to ,7 from 7r + 7r~7 intermediate states for dif- 
ferent cutoff values compared with complete contribution (virtual plus real correc- 
tions, labelled 'inclusive') evaluated in scalar QED (FSR) as well as with contri- 
bution from 7r + 7r~ intermediate state (a) and integrated contribution to aj^ ad ' 7 as 
function of the cutoff E cut (b). 
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is shown in FigEb- Thus a special care has to be taken when imposing cuts 
on the photon energy below 50 MeV, as it might lead to significant shift for 

4. FSR at NLO and tests of the FSR models 

The upgrade of the PHOKHARA event generator to version 3.0 (2^] 
consisted in adding the diagrams from Fig0 where the photon emitted 
from the initial state is assumed to be 'hard', say with energy above 10 
MeV at DA$NE and 100 MeV at B-factories. 




Fig. 7. NLO contributions to the reaction e + e — > it + it 7 from real (both soft and 
hard) FSR emission (AIFSNLO(S+H)) (a) and virtual corrections to the 7r+7r~7* 
vertex (AIFSNLO(V)) (b). 



From the analysis of the corresponding corrections to the e + e~ — * ir + ir~ 
process [33] in the framework of sQED, this contribution is expected to be 
of the order of 1%. However, the emission of the initial photon reduces the 
invariant mass of the 7r + 7r~ (or ir + ir~ r y) system to the p mass with high 
probability due to the peak of the pion form factor at the p mass. As a 
result, this contribution is strongly enhanced in the region of invariant mass 
of the ir + ir~ system below the p resonance, as shown in FigJHJi and FigEK 
for KLOE and B-factory energies respectively. Suitably chosen cuts can be 
applied to suppress these NLO FSR contributions. In FiglHb one can see 
that the standard KLOE cuts ^3E3, which consist of the cuts on pion 
angles, the missing momentum angle and the track mass (M^ r ), keep the 
NLO FSR contribution below 2% with respect to the ISR cross section in 
the whole interesting region of the two-pion invariant mass. Similarly at 
B-factories, applying the track mass cut only for events with Q 2 < m 2 , the 
NLO FSR contribution is kept at a negligible level (FigEb). 

Again, as in the case of LO FSR contributions, the main problem consists 
in the model dependence of FSR. Till now only few tests were performed 
to verify the model for FSR. However, if one aims at a measurement of the 
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Fig. 8. Comparison of the Q 2 differential cross sections for ^/s = 1.02 GeV: 
IFSNLO contains the complete NLO contribution, while 1FSLO has FSR at LO 
only. The pion and photon(s) angles are not restricted in (a). In (b) cuts are 
imposed on the missing momentum direction and the track mass (see text for 
description) . 



0.06 
0.05 
0.04 
J 0.03 
, 0.02 
i 0.01 


-0.01 
-0.02 
-0.03 



e + e — > 7r + 7r 7(7) 



0° < 6,, < 180° 



0° < 6» ± < 180° 



10.52 GeV 



2 a 



5 a 







0.6 0.8 1 

Q 2 (GeV 2 ) 



0.0s 
0.06 
0.04 
11.02 


-0.02 
-0.04 
-0.06 
-0.08 



e + e — > 7r + 7r 7(7) 



1 1 

4s 1 


1 

0.52GeV 






b 


M M < O.HSCcV for Q 2 































■ f 














| i i i > 


* i ? 




























0° < % < 

- 0" < IS-- < 

i i 


180° 

i 













0.6 0.8 1 

QHGeV 2 ) 



Fig. 9. Comparison of the Q 2 differential cross sections for ^fs = 10.52 GeV: 
IFSNLO contains the complete NLO contribution, while IFSLO has FSR at LO 
only. The pion and photon(s) angles are not restricted in (a). In (b) cuts are 
imposed on the track mass for Q 2 < m 2 . 

accuracy below 1% such tests become indispensable. 

In the present KLOE experimental setup, where only four momenta of 
the pions are measured, a possibility to test the hard part of the NLO 
FSR (AIFSNLO(H)) contribution is to look at the dependence of the cross 
section on the missing invariant mass. Completely different effect of the cut 
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on missing invariant mass on ISR (ISRNLO) and FSR at NLO, as shown in 
Fig llUb . provides a powerful tool for testing the hard part of the IFSNLO 
contributions. A measurement of this few percent effect, depending on the 
two-pion invariant mass Q 2 is within reach of the KLOE experiment. 
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Fig. 10. Dependance of the relative IFSNLO contribution on the cut on missing 
invariant mass M 2 (a) and radiative correction factor A extracted from MC data 
(see text for explanation) (b). 

The tests proposed in section for the FSR at LO and the above tests 
of the hard part of the IFSNLO contributions can be further extended to 
the test of the virtual radiative corrections to the FSR of FigEb following 
the procedure outlined below. 

The pion form factor parameters can be fitted to the data using the 
measured radiative return cross section. The cross section can be splited 
into four parts 



d(JRR 

dQ 2 



ISRNLO + FSRLO + AIFSNLO(V + S) + AIFSNLO(H) . (3) 



Having the form factor parameters, one can subtract from the data, rely- 
ing on the Monte Carlo simulation: 1. the ISR contributions, which in- 
volve QED corrections only, 2. the FSRLO contributions tested via charge 
asymmetries, so well known, and 3. the hard part of the IFSNLO correc- 
tions (AIFSNLO(H)), tested as described above, so known with required 
accuracy. The left over virtual plus soft photon corrections to the FSR 
(AIFSNLO(V+S)) can be written in the following way 



AIFSNLO(V + S) 



da 



Born 



dQ 2 



(lnu;-/ + A(Q 2 



(4) 
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where a n stands for the LO e + e~ — > tt + tt~j cross section with photons 
emitted from the initial state only and the In it; is the standard soft photon 
logarithm with a known function / as a coefficient. It is clear that one 
can now extract the function A(Q 2 ) from data and compare it with the 
analytical value of the tested model. The function A(Q 2 ) can be extracted 
separately in the same way for each cut on invariant mass. Moreover, for 
self consistency of the tested model the values of A(Q 2 ) obtained for each 
invariant mass cut should be the same. It provides a nontrivial check of the 
model, as the relative contribution of the AIFSNLO(V+S) part to the total 
cross section does depend on the invariant mass cut. The results of that 
procedure applied to the Monte Carlo data obtained by PHOKHARA 3.0 
are shown in Fig llflb . 

Having all separate ingredients of the cross section one can extract not 
only the pion form factor but also calculate the NLO corrections to the 
muon magnetic moment coming from the ir + ir~j intermediate state. 



5. Conclusions 

Extensive discussion of possible experimental tests of the model depen- 
dence of FSR for the ir + ir~ hadronic final state has been made. The impor- 
tance of the 7r + 7r - 7 final state contribution to the muon anomalous magnetic 
moment has been emphasised and a method has been proposed to extract 
that contributions from the data using radiative return. 
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